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We study the similarities of jet-medium interactions in several quite different astrophysical systems using 2D and 3D 
hydrodynamical numerical simulations, and find many similarities. The systems include cooling flow (CF) clusters of 
galaxies, core collapse supernovae (CCSNe), planetary nebulae (PNe), and common envelope (CE) evolution. The sim- 
ilarities include hot bubbles inflated by jets in a bipolar structure, vortices on the sides of the jets, vortices inside the 
inflated bubbles, fragmentation of bubbles to two and more bubbles, and buoyancy of bubbles. The activity in many cases 
is regulated by a negative feedback mechanism. Namely, higher accretion rate leads to stronger jet activity that in turn 
suppresses the accretion process. After the jets power decreases the accretion resumes, and the cycle restarts. In the case 
of CF in galaxies and clusters of galaxies we also study the accretion process, which is most likely by cold clumps, i.e., 
the cold feedback mechanism. In CF clusters we find that heating of the intra-cluster medium (ICM) is done by mixing hot 
shocked jet gas with the ICM, and not by shocks. Our results strengthen the jet feedback mechanism (JFM) as a common 
process in many astrophysical objects. 
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1 Introduction 



Jets are observed in a large variety of astrophysical objects. 
In many of these objects jets play a significant, and even cru- 
cial, role in the evolution due to their high energy content. 
In some of these systems the effects of the jets are regulated 
by a negative feedback mechanism, the jet-feedback mech- 
anism (JFM). Such is the case in cooling flow (CF) clus- 
ters and CF galaxies, where the jets heat the intra-cluster 
medium (ICM), and during galaxy formation when the jets 
were required to heat and expel the ISM. A feedback mech- 
anism might also operate in the explosion of core collapse 
supernova (CCSN), where it is possible that jets launched 
by the newly formed neutron star (NS) or a black hole (BH) 
explode the star. 

The aim of this report is to emphasize the connections 
and similarities between some of these seemingly very dif- 
ferent kind of objects. These similarities allow us to learn 
from one object on the other. The systems studied in this 
paper are listed in Table [T] The report is concentrated on the 
similar processes of jet-medium interaction in the different 
objects. For that, many works on the different objects and 
processes that do not consider the similarities between the 
different objects will not be mentioned here. 
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1.1 Morphology 

The efficiency of the JFM comes from the deep potential 
well of the jet's launching site. In CF clusters, CF galaxies, 
and during galaxy formation the jets are lunched by a su- 
per massive black hole (SMBH), while in CCSN they are 
launched by the newly formed NS or BH. In all cases the 
size of the compact object is several orders of magnitude 
smaller than that of the surrounding medium, ensuring that 
few percents of the accretion gravitational energy can sub- 
stantially influence the ambient medium. 

However, there are two other conditions for the JFM 
to work efficiently: (1) The jets must deposit their kinetic 
energy in the inner region, where the feedback mechanism 
must work. (2) Radiative loses, usually by photons but by 
neutrinos in CCSNe, must be small. The deposition of the 
kinetic energy is via strong shock waves. The shocked jet's 
material temperature can be much higher than the ambient 
medium temperature and, due to pressure balance, its den- 
sity is lower. If radiative losses are small indeed, then the 
post shock jet's material expands and forms a low-density 
bubble. As usually there are two opposite jets, a structure 
of two opposite bubbles is formed. A pair of bubbles (or 
several pairs) is a generic structure of the JFM when they 
can be observed. This structure is termed bipolar in the case 
of planetary nebulae (PNe). In young stellar objects (YSO) 
there are jets that interact with the dense surrounding gas. 
However, radiative cooling is very rapid for YSO jets, and 
no bubbles are formed. Any JFM, if exist, is of very low 
efficiency. 
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Property / System 


Clusters 


galaxy formation 


CCSNe 


PNe 


CE 


Typical Energy (erg) 


10 60 


10 by 


10 bi 


10 44 


10 44 - 10 48 


Typical Mass (M«) 


10 12 


10 11 


10 


1 


1 


Typical Size 


100 kpc 


10 kpc 


10 9 cm 


0.1 pc 


10 - IOORq 


Typical Time 


10 7 - 10 8 yr 


10 7 - 10 8 yr 


1 -3s 


10 - 100 yr 


1 yr 


^bubble ( K) 


10 9 - 10 10 


10 9 - 10 10 


10 10 


10 6 


10 8 - 10 10 


^"ambiant ( -K) 


10 7 - 10 8 


10 6 - 10 7 


fewxlO 9 


10 4 


10 5 


Compact object 


BH 


BH 


NS or BH 


MS or WD 


NS or WD 


mass (M©) 


10 8 - 10 10 


10 6 - 10 9 


1-5 


1 


1 


Main Effect of jets 


Heating ICM 


Expelling gas 


Exploding the star 


Shaping the PN 


Reducing accretion 


Observations 


X-ray bubbles 


(Massive outflows) 


(Axisymmetry) 


Bipolar PNs 





Table 1 Systems discussed in this paper where feedback and shaping by jets take place. The different listed values are 
typical and to an order or magnitude (or two even) accuracy only. Typical energy: Energy in one jet episode. Typical Mass, 
Size: of the relevant ambient gas. Typical time: the duration of the jets activity episode. In the last row of observations, in 
parenthesis are expected observations. Acronym: PNe: Planetary nebulae; CCSNe: core collapse supernovae; CF: cooling 
flow; ICM: Intra-cluster medium CE: common envelope; BH: black hole; NS: neutron star; MS: Main sequence star 



The most prominent similarity of the bubble-pair (bipo- 
lar) structure is perhaps the very similar morphological 
structures found in s ome PNe and some X-ra y deficient bub- 
bles in CF clusters (iSoker & Biskerl l2006h. Two compar- 
isons a re presented in Figure [T] taken from lSoker & Bisker 
2006) where more examples are given. Despite the sev- 



eral orders of magnitude differences in size, energy, mass, 
and timescales, the similarity is not only in the morphol- 
ogy, but in some basic physical processes as well. These 
similarities were studi ed in a series of papers by one of us 
( Soker & Bisker 2006 and references therein). Point sym- 
metric PNe, like Hb 5 presented in the right panel of Fig- 
ure [T] are thought to be shaped by stellar binary interac- 
tions that caused a jet-precession. The similarity of the mor- 
phology of Hb-5 and the p air of bubbles in t he C F cluster 
MS 0735.6+7421 brought fpizzolato & Sokerl d2005h to sug- 
gest that MS 0735.6+7421 has a massive binary BHs system 
in its center. Simulations of PNe jet- shaping are presented 
in section [2 while jet-inflated bubbles in CF clusters are 
presented in section [3] 



In PNe gravity is negligible and there is no need for a 
feedback mechanism. We note though that in some PNe and 
other nebulae around stars, most of the kinetic energy of the 
expanding gas comes from the jets. This seems to be the 
case in the bipolar neb ula, the Hom unculus, around the bi- 
nary system 77 Carinae dSokeriE004a|) . The relevance of PNe 
is that the bipolar structure of pairs of bubbles is well re- 
solved, and can teach us about the interaction of jets with the 
ambient gas. This is demonstrated in the right panel of Fig- 
ure Q] where a much better resolved image of Hb-5 reveals 
many details that cannot be resolved in MS 0735.6+7421 . In 
some cases the jets and bubbles cannot be observed at all. 
This is the case in jets that might exist in CCSNe and dur- 
ing galaxy formation. Simulations of jet-driven CCSN are 
presented in section [4] 



1.2 Negative Feedback 

There are some common processes and properties of the 
JFM that are common to the different astrophysical objects. 
We list them here, and later implement some of them in the 
numerical simulations. 



1. 



Accretion disk. The jets are launched by an accre- 
tion disk around a compact object. In the case of jets 
launched by SMBH such disks are inferred from obser- 
vations, and the jets are directly observed. However, if 
CCSNe are driven by jets, they cannot be observed. The 
same hold for jets during the galaxy formation epoch, 
as these objects are at large distances and in many cases 
the central region is expected to be obscured. In most 
CCSNe there is no sufficient angular momentum to ex- 
plain a continuous accretion disk. Instead, we assume 
that as a result o f the stationary accretion shoc k insta- 
bility (SASI: e g.lBlondin & Mezzacappall2007l but see 
Nordhaus et al.1 2010l for suppression of this instability 



in 3D simulations), or some other stochastic processes, 
segments of the post-shock accreted gas (inward to the 
stalled shock wave) possess local angular momentum. 
Accretion of dense gas from var ying directions can b e 



seen in the recent simulations by (M ueller et al. . 2012bl) . 



Such an accretion form a jittering pair of opposite jets 
dPapish& Sokerl EoTTl) . 

2. Universal jets' properties. The properties of jets 
launched by the compact objects have some universal 
average properties, (i) The velocity of the pre- shock 
jets' material of vj ~ v esc , where v esc is the escape ve- 
locity from the compact object, (ii) The ratio between 
mass lose rate in the two jets to mass accretion rate is 
77 = M f /M acc ~ 0.1. 

3. Non-penetrating jets. For a jet to deposit energy in the 
relevant inner regions it should not penetrate through 
these regions; this is termed the a non-penetrating JFM. 
The condition for jet- stopping is that the time required 
for the jet to propagate through the surrounding gas 
and break out of it must be longer than the typical time 
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Perseus 

(Cluster of galaxies; 
X-ray) 

(Fabian et al. 2000) 



NGC 3587 

(Planetary nebula; visible) 



(Guerrero et al. 2003) 




MS0735.6+7421 
(Galaxy cluster; 
X-ray + radio contours) 

(McNamara et al. 2005) 



Hb5 

(Planetary nebula; visible) 



(Terzian & Hajian 2000; Schwarz et al. 1992) 




Fig. 1 Comparing false color X-ray image of a galaxy 
cluster with a visible image of a planetary nebula (PN). 
The upper panel emphasizes pairs of fat spherical bub- 
bles near the center, while the bottom panel empha- 
sizes point symmetric morphology of the bubble pair, 
and the consequence s of bad re s olution (right PN im- 
age). I mages from [Fa bian et al. d2000h. Guerrero et al. 
d2003k. iMcNamara et alj d2005h . iTerzian & Haiianl d200< ' 
and lSchwarz et al.1 dl992l) . 

www.an-journal.org 



for restarting the penetration process. The restarting 
can be due to transverse motion of the jets, as in the 
jittering-jet model or in the case of jet precession, due 
to the transverse relative motion of the surrounding 
gas and the SMBH (hence the jets continuously en- 
counter fresh gas), due to the orbital motion of the 
companion in a common envelope (CE), and by wide 
jets that propagate very slowly through the ambient 
gas. These cases have been explored in t he pa st for jets 



launched by SMBHs. ISternberg etaTI d2007k showed 



that slow massive wide (SMW) jets can inflate the fat 
bubbles that are observed in many CFs, in clusters, 
groups of galaxies, and in elliptical galaxies. The same 
basic physics that prevents wide jets from penetrating 
through the ICM was sho w n to h ol d for precessing jets 
(Fal ceta-Goncalves et all 1201 Ol: ISternberg & Sokej. 
2008a|) . or a relative m o tion o f the jets to the 



medium dBriiggen et all 2007 : Mendvgral et all 



120121: iMorsony et all l2010l: Isokek [2009b. If the"jets 



penetrate to a too large distance, then no bubbles 
are formed, while in intermediate cases elongated 
and/or detached from the center bubbles are formed 



2003 



2012; 



Alouani Bibi et al. 2007; Basson & Alexander 


Heinz et al. 


l2006t iMendvgral et al. 


201 ll 

20 id 


lOmma et al 


J 12004]: O'Neill & Jonesj 



Sternberg et al. 2007; Vernaleo & Reynolds 2006b . 



4. Available mass. The mass available for the inflow is very 
large. Namely, the mass that is eventually accreted to the 
NS in CCSN or the SMBH in CFs is limited by the JFM 
and not by the mass available in the surroundings. 

5. Stochastic accretion. In the JFM considered here in 
CCSN and CFs (both in clusters and in galaxies), the 
feeding of the accretion disk is by cold clumps. By cold 
we refer to a gas temperature much below the virial 
one. Although the average rate of accretion is regulated 
by the feedback process, there might be large tempo- 
rary variations in the rates of mass and angular momen- 
tum accretion values. The cold accretion allows both 
accretion and jet activity to coexist. As well, the cold 
clumps falls quite rapidly, allowing fast communication 
between the ambient gas and the accretion disk, as re- 
quired in the JFM. In CFs this is termed the cold feed- 
back mechanism, and it is discussed in section 13.11 In 
CCSN the stochastic accretion of clumps allows a tem- 
porary formation of accretion disk. 

6. Incomplete energy deposition. The JFM we explore here 
differs from other mechanisms by the requirement that 
its influence on the environment is not complete. In CF 
clusters and galaxies, the accretion is of cold clumps; 
the cold feedback mechanism. Namely, the heating by 
the jets although very efficient, is not complete. Most of 
the cold gas forms stars and cold filaments in the ICM, 
rather than feeding the SMBH. The small fraction of the 
cold gas mass that does feed the SMBH is sufficient to 
maintain the JFM. Hence, star formation and cold fila- 
ments in CF clusters and galaxies are generic features of 
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the JFM considered here. In CCSN the accretion contin- 
ues for a few seconds while the jets and the bubbles they 
inflate expel most of the core mass. The regulation of 
this continuous accretion episode is part of the feedback 
mechanism. This is in contrast to most models of CC- 
SNe explosion that start by reviving the accretion shock 
and terminate accretion, and where a feedback mecha- 
nism does not exist. 

2 Bipolar shaping of planetary nebulae 

There are a number of numerical simulations that try 
to reproduce the asymme trical structure of PNe (e.g., 
Huarte-Espinosa et al. 1201 2l for a recent paper and refer- 
ences therein). We here concentrate only on those PNs that 
have two pairs of bubbles similar in structure to pairs of bub- 
bles in CF clusters. These are most likely formed by jets. 

Our simulations are performed by usi ng version 4.0-beta 
of the FLASH code dFrvxell et alll2000b . The FLASH code 
is an adaptive-mesh refinement modular code used for solv- 
ing hydrodynamics or magnetohydrodynamics problems. 
Here we use the unsplit PPM (piecewise-parabolic method) 
solver of FLASH to simulate gas dynamics in different 
astrophysical environments. The simulations are differ by 
adding different cooling functions and by the appropriate 
EOS. In the simulations of PN shaping we include radia- 
tive cooling, but not gravity. In the simulation of CCSN de- 
scribed later we do incorporate gravity. 

The radiative cooling is added to the simulation at all 
temperatures T > 10 4 K, and it is carefully treated near 
contact discontinuities to prevent large temperature gradi- 
ents from causing unphysical results. The cooling func- 
tion for solar abundances that we use was taken from 
dSutherland & Dopita lfl993l: their table 6). 

We employ a full 3D adaptive mesh refinement (AMR, 
6 levels (2 9 cells in each direction)) using a Cartesian grid 
(x,y,z) with outflow boundary condition at all boundary 
surfaces. We define the x — y (z=0) plane to be the equatorial 
plane of the PN and simulate the whole space (the two sides 
of the equatorial plane). The grid size is 2 x 10 17 cm in the 
x and y directions, and 4 x 10 17 cm in the z direction. 

At the beginning of the simulation, t = 0, the grid is 
filled with a spherically- symmetric slow wind, blown by 
the AGB stellar progenitor, having a uniform radial veloc- 
ity of v W i n d = 10 km s _1 . The density at t = is taken 

where we take here a slow- 
= l(r 5 M yr" 1 . The jet 



to be p(t = 0) - 

wind mass loss rate of M W i n d 
is lunched from the inner 6 x 10 15 cm region, and within a 
half opening angle of a = 70 °. This is a wide jet; similar re- 
sults will be obtained by rapidly precessing jets. The jets are 
launched with a radial velocity of Vj et — 600 km s _1 and 
the mass loss rate in each side of Mj e t = 3 x 1O _8 M yr _1 . 
For numerical reasons a weak slow wind is injected in the 
sector a < 6 < 90 °. The slow wind and the ejected jet 
start with a temperature of 10000 K. The initial jets' tem- 
perature has no influence on the results (as long it is highly 



supersonic) because the jet rapidly cools due to adiabatic 
expansion. 

Our calculations do not include the ionizing radiation 
and the fast wind blown by the central star during the PN 
phase. We simply aim at showing the shaping of the nebu- 
lar gas to a bipolar shape by blowing wide jets. There is no 
feedback heating, but there is a shaping that is very similar 
to the one in clusters. The mass of the hot bubble is very 
small relative to the ambient mass, while the linear momen- 
tum of the jet's material is comparable to the one of the AGB 
wind. 

In Figure [2 we show color-maps of the the gas density, 
enlargement of the vortex zone, and the gas temperature, 
at t = 580yr and t = 1140yr, in the upper row (panels 
[(a)| , |(b)[ and[(c)|), and lower row (panels |(d)| , [(e)! and [(f)]), 
respectively. In all plots we show the meridional plane of the 
nebula. We recall that the simulation is full 3D, and include 
both side (no mirror symmetry is assumed in this section). 

The following features should be noticed. 

1. The two opposite jets form a bipolar structure of two 
'fat bubbles'. By fat bubbles we refer to more or less 
spherical bubbles attached to their origin. We will find 
this structure in the simulations of bubbles in clusters of 
galaxies in the next section. 

2. A dense thin shell is formed around the bubbles. This is 
also seen in bubbles in clusters of galaxies. 

3. A structure of two shocks and a contact discontinuity 
between them is formed. The forward shock runs at a 
very low Mach number into the slow AGB wind. The 
jet is shocked in the reverse shock to very high temper- 
atures. The result is that the bubble is much hotter than 
the shocked ambient gas. The same holds in bubbles in 
clusters, but not in bubbles in CCSNe (see Table [T]). 

4. Vortices are formed on the sides of the bubbles. We will 
encounter much prominent vortices in our 2D simula- 
tions of bubbles in clusters. 

5. The vortices to the sides of the bubbles compress gas 
to a thin disk in the equatorial plane. This is best seen 
in panels |(b)| and [(e)J The vortices not only form the 
disk-like structure, but also slow it do wn. This process 
was p reviously studied numerical by lAkashi & Soked 
d2008h . who also suggested that part of the equatorial 
dense gas might flow back to the center and form a large 
(up to hundreds of AU) Kepleria n disk; Planet s might be 
formed in such Keplerian disks dPeretsl.l2Q10h . 

3 Jet feedback in galaxies and clusters 
3.1 Accretion 

It is widely accepted that feedback powered by active galac- 
tic nuclei (AGN) has a key role in galaxy formation and 
in cooling flows (CFs) in galaxies and in clusters of galax- 
ies. In gal axy formation AGN feedback heats and expels 
gas (e.g.. lOstriker et al.ll2010l and references therein), and 
by that can determine the correlation between the central 
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Fig. 2 Color-maps of the gas density (in g cm -3 ), density and velocity in an enlarged vortex zone, and the gas temper- 
ature (in K). Upper and lower rows are for t = 580 yr and t = 1140 yr, respectively. The velocity arrows in panels [(a)! 
and |(d)| show only the direction, while in the panels with enlarged vortex zones they are linear with the magnitude of the 
velocity. 



SMBH mass and some properties of the galaxy. In cooling 
flow clusters jets launched by the SMBH heat the gas and 
maintain a small, but non zero cooling flo w (see review by 
Fabianll2012l: iMcNamara & Nulsenll2012h: this is termed a 



moderate cooling flow in Sok er et al . (2001). 



There is a dispute on how the accretion onto the 
SMBH occurs, in particular in cooling flows. One side 
argues for accretion to be of hot gas via the Bondi ac- 
cretion process (e.g 
l201ll: iRussell et al. 



Allen et all 120061: iNaravan & Fabiani 
201(1 



while the other side argues 



that the accretion is of dense and cold clumps in what is 
termed the cold feedback mechanism (IPizzolato & Soker . 



20Q5|) . The cold feedback mechanism has been strength- 



ened recently b y observations of cold gas and by more de_ 
tailed studies dCavagnolo et al. ^ 20 1 1 ^ Farage et all 20 1 2 : 
[ Gaspari et all l2012allb|: iKashi et all |2012l: iMcCourt et al. ' 
20121: iNesvadba et all BOllt IPop« " 



2009; Revazetal 



20081: ISharma et all 2012: Wilm an et all 120091 1201 lh 



The Bondi accretion process , on the other hand , 
suffers from several problems (|Cavagnolo et all 12011 
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1 ° „15 



pc] 



Fig. 3 The average density, pressure and temperature of 
the hot bubble as a function of radius. There are two free 
parameters: the fraction of the escape velocity in which the 
hot gas leaves the hot bubble, and the fraction of mass from 
the winds that is incorp orated into the bubble . The parame- 
ters used here are from (Hillel & SokeilEoi^ . The pressure 
of the hot bubble at r — 20 pc about equals the ISM pres- 
sure of NGC 3115. 



McNamara et all 1201 ll: ISokeii l2009h . We point out yet 
another problema tic point with it. In a recent paper, 
Wong et al.l (1201 1|) resolved the region within the Bondi ac- 
cretion radius of the SO galaxy NGC 3 1 15. If the density and 
temperature profile is interpreted as resulting from a Bondi 
accretion flow into the M B h = 2 x 1O 9 M central SMBH, 
the derived accretion rate is Mb = 2.2 x 1O _2 M yr _1 . 
They note that for a radiation power of O.lM^c 2 , the 
expected accretion luminosity is six orders of magnitude 
above the observed upper limit. They attribute this to a pro- 
cess where most of the inflowing gas is blown away. 

We take a different view. We argue that the Bondi accre- 
tion flow is not relevant for the conditions in typical galaxies 
and clusters of galaxies. The reason is that one cannot as- 
sume a zero pressure at the center, either because of stellar 
winds or because of jets and winds blown by the AGN. 

In application t o the galaxy NGC 3115, we 
(IHillel & SokerL |2012|) have estimated the average density 
and pressure of a hot bubble around the SMBH. The hot 
bubble is formed by stellar winds. In Fig. [3 the average 
density, pressure and temperature of the hot bubble are 
shown as a function of radius. The pressure of the shocked 
stellar winds of the high- velocity circum-SMBH stars can 
be larger than the ISM pressure near the center. This ac- 
counts, we argue, for the accretion rate of NGC 3115 being 



much lower than the Bondi accretion rate (I Wong et al 
2011b . 



The average density of the hot shocked stellar wind 
within the bubble of radius R is given by 



and the pressure is 

2 /4tt 
3~ 

R 



R c 



R 



Pu*{R) 



I 



-n*(r)r]rn*ul (r) 



(2) 



47rr 2 dr. 



Here, u* (r) is the velocity of a star at a distance r from the 
SMBH, (3u*(R) is the velocity at which the hot gas leaves 
the hot bubble, m* is the average mass loss rate per star, 
n* (r) is the number density of stars near the SMBH, and 
r] is the fraction of mass from the stellar winds that is in- 
corporated into the bubble. Fo r the stellar number de nsity in 
NGC 31 15 we take (based on lKormendv et al.lll996h 



,^ _s J 1, r < 3 pc 

n* = 5 x 10 5 pc 3 J 



[(r/3pc)- 3 , r > 3pc. 



(3) 



In conclusion, we find that in some cases this pressure 
exerted by the stellar winds is significant and can substan- 
tially suppress the inflow of the ISM relative to what a sim- 
ple Bondi accretion would give. Shocked winds of circum- 
SMBH high- velocity stars form a bubble of hot gas whose 
pressure is significant. There are some uncertainties in the 
model, such as the exact behavior of the stellar mass loss, 
trajectories of stars around the SMBH, and the stochastic 
behavior of the post-shock stellar winds. Some of these will 
be studied in future numerical simulations. However, the re- 
sult that the stellar winds cannot be ignored is robust. 

For some values of the parameters we found that a situa- 
tion might arise where the hot bubble's density is lower than 
the ISM density. In this case Rayleigh-Taylor (RT) insta- 
bility takes place, and a density-inversion layer is formed. 
Such a structure, we claim, is similar to the density inver- 
sion found in the outer atmosphere of r ed giant stars (e.g., 
iFrevtag & Hofned 120081 lHarpazI fl984h . but not identical. 
Although hot tenuous gas buoys outward and dense ISM 
gas moves inward, the density-inversion layer itself contin- 
ues to exist. The ISM gas is heated near the center and ac- 
cumulated to the hot bubble. 

Our result is more general in showing that in many cases 
the Bondi accretion process does not work because one 
of its basic assumptions, that there is no central pressure, 
breaks down. This is one of several reasons why the Bondi 
accretion model is not applicable in many cases. 

We note that the considerations above cannot be applied 
straightforward to galaxies like the Milky Way where the 
SMBH mass is relatively low. The reason is that the mass 
of SMBH at the center of the Galaxy is much smaller than 
that in NGC 3115, while the stellar velocity dispersions a 
are similar, as well as the stellar densities. The radius of 
influence of the SMBH, GM Bil /a 2 , in NGC 3115 is large, 
~ 30 pc, and contains many stars with significant total mass 
loss rate that can act against the ISM pressure. In the Milky 
Way the radius of influence is only ~ 0.1 pc. There are 
many stars with large mass loss rate outside this regions. 
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There is continuous mass loss rate by stars from large dis- 
tances inward. Thi s is seen in the numerical simulations of 
Cuadra et al. 2006. We note, though, that when they take a 
smaller inner radius in their simulations, the accretion rate 
decreases. We interpret this as the building of pressure near 
the center. 

3.2 Bubbles inflation 

Bubbles (cavities) devoid of X-ray emission, mostly as op- 
posite pairs, are observed in a large fraction of cooling flow 
(CF) clusters and gro ups of galaxies, a s well as in CF el- 
liptical galaxies (e.g jDong et al.lEoio|) . These bubbles are 
inflated by jets launched from the central active galactic 
nuclei (AGN), as evident by the radio emission that fills 
most bubbles. In some cases two opposite chains of bubbles 
that are close to each other, and ev en overlap, are observed, 



as in Hydra A (I Wise et all |2007|). and in two bu bbles in 



the galaxy group NGC 5813 dRandall et all 1201 ll) . These 
chains were usually attributed to several episodes of jet ac- 
tivity. 



Refaelovich & Sokeri (12012 ) suggest another plausible 



mechanism for cavity chains creation. Several fragmenta- 
tion mechanisms act on the primary vortex that created 
just behind the jet's head and split it to several smaller 
vortices. These mechanisms include Kelvin-Helmholtz and 
Rayleigh-Taylor instabilities which the jet's cocoon flow is 
prone t o. Using hydrodynamical sim ulations with a contin- 
uous jet lRefaelovich & Sokeri (|2012|) show that such a frag- 
mentation process takes place for a wide range of jet pa- 
rameters, and causes the fragments to appear like a chain of 
cavities in the X-ray brightness map. 

In the current paper we present new results which may 
help to discriminate observationally between the two mech- 
anisms, vortex fragmentation and multiple jet-launching 
episodes. 

To emphasize the differences between the above men- 
tioned me chanisms we repeat the stan dard simulation pre- 
sented by Refa elovich & Soker (2012) with the difference 
that the jet is switched on and off. The activity phase last 
20 Myr followed by a pause of 50 Myr. We simulate four 
such cycles spanning a period of 280 Myr. The total two 
jets power is P^j = 4 x 10 45 erg s -1 , the jets Mach num- 
ber is 15 relative to the ambient medium, their velocity is 
Vj = 1.33 x 10 4 km s -1 , and the half-opening angle of 
each jet is a = 30°. 

We use the PLUTO hydrodynamic code 
(Mign one et all l2007|) . The simulations are 2.5D in 



the sense that we use a spherical coordinate system, but 
impose cylindrical symmetry, hence calculating the flow 
with a 2D polar grid (r, 6). 

In Figures 01 and \5\ we present the results of the two 
simulations mentioned above. One run with the continuous 
jet and another with the jet switched on and off. Figure 
shows density color-maps and Figure \5\ shows their simu- 
lated X-ray map counterparts after applying an unsharp fil- 



ter as described in Refaelovic h & Soker (2012). The figures 
compare the two runs on three different stages. 

1. t = 20 Myr - Just before the end of the first ac- 
tivity episode in the switched run. The positions of 
the shocks are identical, the exact cavity structure 
may be somewhat differen t due to convergence issue 
dRefaelovich & Sokeri 120121) . 

2. t = 69 Myr - The end of the first pause episode, i.e., 
just before the next activity starts. 

3. t = 180 Myr in the switched run Vs. t = 100 Myr in 
the continuous run. In the switched run it is 20 Myr after 
the end of the 3rd jet launching episode. 

There are several conclusions following these simula- 
tions. 

- Cavity position relative to the shock front. During the ac- 
tive phase the primary bubble (cavity) which is created 
near the head of the jet is located just behind the shock. 
The jets' momentum transferred to the bubble keeps the 
bubble behind the shock. During the idle phase, the jet 
stops driving the bubble and the shock separates from 
the cavity and runs ahead (figs. |4(c)| and |5(c)| ). This em- 
phasizes the fact that the outward bubbles' motion is 
driven mainly by the momentum of the jet, rather than 
by buoyancy. It should be noted that the same process of 
jet driving and separation will occur with each next jet 
activity episode. 

- Positions of secondary shocks. Secondary shocks from 
successive activity episodes catch-up with the previous 
shocks very fast. However, they catch-up only in the 
direction of the jet. In the perpendicular direction the 
shocks stay separated for much longer times (they actu- 
ally do not catch up during the entire simulation time). 
This point is illustrated in Figure |4(e)| and its X-ray 
counterpart in Figure [5(e)) There are two shocks which 
catch-up at z = 250 kpc on the symmetry axis, fol- 
lowed by a third shock around z = 125 kpc. It may be 
seen that all the three shocks are well separated along 
the R direction (note that the outermost shock is par- 
tially outside the frame). 

Thus, while it is probably impossible to count the 
number of episodes using the shocks on the chain 
axis, it should be possible, at least in principle, to 
discr iminate between the m at the perpendicular direc- 



tion. 



Ran dall et al.1 d201l|) find three such shocks in the 



galaxy group NGC 5813, and argue for multiple jet- 
launching episodes to account for the chain of bubbles 
there. 

- Inter-cavity separation. When considering the separa- 
tion between X-ray deficient cavities we should con- 
sider two different cases. In the first case the recent 
shock didn't catch up yet with the bubble inflated by 
the previous activity episode. At this stage we expect to 
see very sharp separation between the cavities because 
of the shock front that runs between them as in Fig- 
ure [6] Note that the outermost bubble is detached from 
the shock and therefore not sharply separated from the 
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Fig. 4 Colormaps of electron density in logarithmic scale and units of cm -3 . Results of the continuous ( |(b)|(d)| 
and switched ( |(a)|(c)[[(e)| ) runs at different times. Axes are in kpc. The horizontal boundary (labelled Z) is the symmetry 
axis of the two jets. The left vertical boundary of the figure (labelled R) is taken along the mirror- symmetry plane of the 
flow z = 0; it is termed the equatorial plane. The jet is injected at r = 1 kpc, a region that is not well resolved in the 
figure. In panel |(b)| there are velocity arrows depicting the vortices present inside the bubbles. There are 3 velocity bins: 
500 - 1000 km s~\ 1000 - 5000 km s _1 and 5000 - 10 4 km s _1 - from shortest to longest. Velocities below 500 km s _1 
and above 10 4 km s _1 do not appear in the figure for clarity. 
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Fig. 5 Synthetic X-ray image of the simulations viewed from the equatorial plane direction (inclination of 90°) per- 
formed by integrating n 2 e along the line of sight. Results of the continuous ( |(b)|(d)| and switched ( |(a)|(c)| [(e)| ) runs at 
different times. 



ICM in the X-ray map. The second case occurs when a 
newer shock interacts and disrupts a previously inflated 
bubble. At this case the cavity (which consists of a vor- 
tex) is fragmented to smaller vortices which may appear 
as additional cavities. No sharp separation between the 
cavities exists in this case, and the image will resem- 
ble the multip le cavities formed by a con tinuous jet as 
described by iRefaelovich & Sokerl (I2012|) . We want to 
emphasize that both cases can in principle coexist in the 
same system. For example, in a case of three episodes 
when the second shock already disrupted the first cavity 



but the third shock didn't reach yet the second cavity as 
in Figures [4(e)| and [5(e)) 

3.3 Heating by mixing 

To study the heating of the ICM by jets we performed 2.5D 
simulations with the code PLUTO. We launch jets with a 
wide half-opening angle of 6 — 70° and a velocity of 
9600 km s _1 . To follow the lasting effects of the jets we 
run the simulations for a significantly longer time beyond 
the jet activity phase. 
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Fig. 6 Electron density color-map and its X-ray map counterpart from the switched run (as in Figures |4] and 0). This 
figure shows the run 3 Myr after the second jet activity ceased. The first shock is at z = 175 kpc and the second shock 
is at z ~ 75 kpc. The two bubbles are clearly seen. This figure depicts the strong cavity separation which occurs when 
the shock from a secondary jet activity episode is located between the newly inflated bubble and its predecessor from a 
previous activity episode. 




t=isoMyr 



Figure[7]shows the evolution of the density in the merid- 
ional plane {vd,z) at three times and for three models. 
The left and middle columns show results for a single jet- 
launching episode lasting 20 Myr, and with a kinetic two- 
jets power of 2 x 10 44 erg s _1 and 2 x 10 45 erg s _1 , respec- 
tively. The right column shows results for two subsequent 
jet-launching episodes, both lasting 20 Myr and with a ki- 
netic two-jets power of 2 x 10 45 erg s _1 ; the first episode 
starts at t = and ceases at t = 20 Myr, while the sec- 
ond one starts at t = 60 Myr and ceases at 80 Myr. In each 
panel the z (vertical) axis is the symmetry axis (the jet axis). 
Only one half of the volume is simulated. 

The first row shows the flow structure at the end of 
the (first) jet injection episode. The inflated bubbles are 
clearly seen. At early times they are attached to the cen- 
ter and form 'fat bubbles'. It is evident that a lower power 
jet (left panel) inflates a smaller bubble. By t = 85 Myr 
in the left and middle columns the low-density bubble has 
risen, leaving behind a trail of vortices. In the right col- 
umn t = 85 Myr corresponds to 5 Myr after the second 
jet episode has ceased. The shock front excited by this sec- 
ond jet-launching episode is not spherical, contrary to the 
case with the first jet episode. The third row shows the flow 
structure at the end of our simulations, long after the end of 
any jet activity. 

In Figure [8] we plot velocity arrows that emphasize the 
turbulent nature of the bubbles and their surroundings. Vor- 
tices are clearly seen, and can be compared to those studied 
in section lT2l This turbulent flow is responsible for the mix- 
ing of very hot shocked jet material with the ICM. 

Figure [9] shows the synthetic projected X-ray maps (ob- 
tained by the same procedure as described in section IX2l) of 
the high-power runs at t = 85 Myr (corresponding to the 
middle row of the middle and right columns of Figure [7]). 
We mirrored the space simulated twice (about the symmetry 
axis and about the equatorial plane) to obtain the full image. 
One pair of bubbles away from the center is clearly seen 




25.5 



25 50 
td [kpc] 



26.5 



Fig. 8 The density map in the (vj,z) plane for the simu- 
lation with one jet episode and a one-jet kinetic power of 
E = 10 45 ergs -1 . Also shown are velocity vectors, for 
v > 120 km s -1 . 



in the left panel. The right panel shows two pairs of bub- 
bles that were generated by the two jet-launching episodes. 
These synthetic maps can be compared with those in Figure 
[5] that were obtained for different physical parameters, e.g., 
narrower jets. Rich variety of bubble morphologies (X-ray 
deficient bubbles) similar to those observed can be formed 
with different parameters. 

Figure [TO] shows the evolution of the average temper- 
ature of gas (the "TR91 tracer gas") that was initially in 
a circular region in the (w,z) plane (torus in 3D), cen- 
tered around (m, z) = (9, 1) kpc, and having a radius of 
0.25 kpc. We follow this gas using an artificial flow variable 
£, which is frozen-in to the flow. We denot e this gas TR91 
(for more details see lGilkis & Sokern2012|) . Four cases are 
shown, for two values of jet power, and for one or two jet 
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Fig. 7 The density maps in the meridional (w, z) plane at three times (as marked on each panel), and for three cases 
(as marked above each column). The first number is the power of one jet in erg s _1 , followed by the active jet period(s) 
in Myr: given inside each square parenthesis is the time the jet starts and ends. The density scale is logarithmic, from 
20-26.5 g cm -3 in blue to 10 -24 5 g cm -3 in red. The symmetry axis of the jets is along the vertical (z) axis. 
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Fig. 9 The projected X-ray map in the full (vj,z) plane 
at t = 85 Myr for two cases (as marked above each panel). 
The first number is the jet power in erg s~ 1 , followed by the 
active jet period(s) in Myr. The X-ray image is obtained by 
integrating over the density squared. 



episodes, as marked in the inset of the figure. Clearly seen 
are sharp rises in temperature that occur when the forward 
shocks hit the gas in this tracer (marked on the figure). Fol- 
lowing these shocks are adiabatic oscillations of the temper- 
ature caused by sound waves. Eventually, the temperature 



returns to a value very close to its starting value; the heat- 
ing by the shock has no lasting effect. The only significant 
and lasting heating happens when the tracer material mixes 
with hot low-density shocked jet-material. This happens for 
the cases of the high-power jet(s) (Eij = 10 45 erg s _1 ) at 
t ~ 60 — 90 Myr. The mixing of the TR91 gas is shown in 
Figure fTTl 

The main conclusion from these simulations is that the 
turbulent motion excited by the inflation process of the bub- 
bles mixes the ICM with the shocked jets' material, leading 
to the heating of the ICM. In the next subsection we further 
show that this process is more important than shock waves 
in heating the ICM. 

3.4 Shock heating 

Consider a short duration explosion, namely, a short- 
duration bubble-inflation episode, with a total energy in- 
side a constant density po of the ICM. The constant energy 
medium holds only in the inner few kpc, but the assumption 
is adequate enough for the present goal. The shock radius 
and velocity for an ICM with 7 = 5/3 reads (e.g.. iTayloi] 
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Fig. 10 Average temperature of the TR91 tracer gas for 
four cases as marked in the inset: the first number is the 
jet power in erg s _1 , followed by the active jet period(s) 
in Myr. The TR91 tracer gas is the gas that was located 
initially at z) = (9, 1) kpc. N ote that the low-powe r 
jets do not manage to heat TR91 (see lGilkis & Sokej2012h . 
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Fig. 11 The concentration of jet material £ je t (fraction 
of jet material at each point) is shown by the color coded 
map in logarithmic scale, for the simulation with one jet 
episode and a one-jet kinetic power of E = 10 45 erg s _1 . It 
is clearly seen that the jet material mixes with the ICM. We 
also follow the material of two tracers, TR61 and TR91. 
Here TR91 is the gas that started inside the circular re- 
gion centered on (w, z) = (9, 1) kpc and having a radius 
of 0.25 kpc, with a similar definition for TR61. The black 
(white) contours show where the concentration of TR61 
(TR91) is one percent. Both tracers are well mixed with 
their surroundings. 



R s (t) = 1.15 



0.46 



Pot 3 



1/5 



1.16 



Po 



1/5 



nl/2 



M(R S ) 



(4a) 
(4b) 



where M(R S ) = poV is the ICM mass inside radius R s , 
and V = 4irRg/3 is the volume. The shock Mach number 
M s = where C s = (jP/po) 1 / 2 is the sound speed of 
the ICM, can be evaluated as follows. For the constant den- 
sity and temperature assumed here the pressure is constant 
as well. This is of course not accurate, but we can think of a 
mass-weighted averaged pressure. The total thermal energy 
of the gas within radius r is .E t h = (3/2)PV, and the Mach 
number can be written as 



M, = 1.1 



1/2 



The energy equation reads 

at e t h 



^cool 



(5) 



(6) 



where e t h is the thermal energy per unit volume, L is the 
cooling rate per unit volume, and t coo \ is defined here as 
the cooling time in constant volume (at a constant pressure 
the cooling time is longer by a factor of 5/3). For a change 
A s = A In Pp _5//3 in each shock, the number of shocks re- 
quired within a cooling time to maintain the gas at the same 
entropy is N s (t coo \) = A7 1 . The total energy required by 
the AGN activity during the average cooling time within ra- 
dius r is given by 



F (t \ M * Eth 
E AGN (t cool ) = — — . 



(7) 



This gives E AGN (t coo i) = 112E th for M s = 1.3, and 
E ag N(t cool) = 27£ t h for M s = 1.7. This is a very in- 
efficient heating mechanism. The power of bubbles in clus- 
ters of galaxies is in the order of mag nitude of what is re- 
quired to maintain the hot ICM (e.g., Ma et al.ll2012l for a 
recent reference), namely, ^agat(^cooi) ~ Eth, and cannot 
account for heating by shocks. 



In the g roup of galaxies NGC 5813 iRandall et al. 



(1201 lL |2012|) detected 3 shocks perpendicular to the axis 
of the three bubble pairs. The shocks' parameters are 
[R S ,M S , A-\t cool (R s )} = (1.2 kpc, 1.71, 10, 2 x 10 8 yr), 
(10 kpc, 1.52, 20, 9 x 10 8 yr), and (25 kpc, 1.3, 77, 2 x 
10 9 yr). Namely, to heat the gas at 10 kpc over a time of 
~ 10 9 yr the AGN energy release in jet activity should be 
~ 20 times the total energy within radius of 10 kpc, or 77 
times within a radius of 25 kpc. This is problematic accord- 
ing to our findings for the following reasons. (1) The heat- 
ing along the jets direction is very efficient as mixing takes 
place there very efficiently. This implies that the activity that 
required to keep the gas residing perpendicular to the jets' 
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axis hot by shock heating, will lead to very high tempera- 
tures along the jets direction. (2) The shock-heating feed- 
back requires many bubble-inflation episodes. The bubbles 
remove gas from the center. This will cause gas to flow in- 
ward, hence being more prone to heating by mixing. One 
cannot treat the gas as sitting in one place. A huge circular- 
ization of gas will star. 

Our preferred explanation is that (a) mixing occurs per- 
pendicular to the jets axis close to the center, (b) At larger 
distances the cooling time is longer. Over such a time the 
bubbles directi on is very likely to change, e.g., the Perseus 
ghost bubbles dFabian et all l2000|) . That varying jets' axis 
lead to a much more efficient interaction with the ambient 
gas is discussed in the next section where the explosion of 
CCSN is discussed. 

4 Exploding core collapse supernovae 

(CCSN) 

In this section we describe numerical results for the 
jittering-jet model for CCSNe. 



4.1 Jet-driven explosion 

Jet-driven superno va explosion mod e ls hav e a long his 
tory (e .g. iBisno vatvi-Ko gan et al 
20091 1201 ll iHoflich et alJ 12001 



1999; 



1970; 



Lazzati et al 



Wooslev & Janka 



MacFadven et al 

J" 



2011 
12001 



97(1 

iKhokhlov et al.l 



Couch et al.l 



explosion (Pap ish & Soke 




iLeBlanc & Wilsonl 
iMeier et al.l 1 19761: 

j itteri ng-jet model for 



2012) is based on the 



following points, that differ in several ingr edients from the 
mode ls cited above (for more detail see Papish & Soker 



2011). 



(1) We don't try to revive the stalled shock. To the contrary, 
our model requires the material near the stalled- shock to 
fall inward and form an accretion disk around the newly 
born NS or black hole (BH). In the present simulations 
we do not include the stalled shock and don't resolve the 
accretion disk, but we do show the accretion inflow through 
the equatorial plane. 

(2) We conjecture that due to stochastic processes and 
the stationary accr e tion shock instability (SAS I; e.g. 



Blo ndin et all 120031: iBlondin & Mezzacappal 12007b seg- 



ments of the post-shock accreted gas (inward to the 
stalled shock wave) possess local angular momentum. 
Fogl izzo et al.l ( 2012 ) found in their lab experiment that 
although the inflow contains no net angular momentum, the 
non-linear interactions in the inflow ultimately favor a sin- 
gle spiral mode to be formed. This leads the accreted mass 
to temporarily posses angular momentum. This accretion of 
angular momentum is crucial to our model, as it probably 
leads to the formation of an accretion disk with rapidly 
varying axis direction. We note that in order to revive the 
stalled shock, as is required in neutrino-d riven models, 
the SASI seems to be very important (e.g., iMueller et al 



2012ah . The SASI 
simulations (e.g. 



is 



obtained in n ew high-resolution 
20123) . ms to be 



Mueller et al 



ubiquitous in CCSNe. 

(3) We assume that the accretion disk launches two opposite 
jets. Due to the rapid change in the disk's axis, the jets can 
be intermittent and their direction rapidly varying. These 
are termed jitter ing jets. 



(4) We show in iPapish & Sokeri (120111) that the jets pene- 
trate the inf ailing gas up to a distance of few x 1000 km, 
i.e., beyond the stalled- shock. Beyond fewxlOOOkm 
the jittering jets cannot penetrate the gas any more. The 
jittering jets don't have the time to drill a hole through 
the ambient gas before their direction changes; they are 
shocked before penetrating through the ambient gas. This 
condition can be met if the jets' axis rapidly changes its 
direction. This process of depositing jets' energy into the 
ambient medium to prevent further accretion is termed 
the penetrating jet feedback mechanism. The jets deposit 
their energy inside the star via shock waves, and form hot 
bubbles. These bubbles and their similarity to hot bubbles 
in other astrophysical objects is the focus of this section. 

(5) The jets are launched only in the last phase of accretion 
onto the NS. For the required energy the jets must be 
launched from the very inner region of the accretion disk. 
In our jittering-jets explosion model the jets are launched 
close to the NS where the gas is neutron-rich (e.g., 



Kohri et al. 112005b . As long as the jets and the bubbles they 
form don't explode the core, accretion continues. When the 
jets manage to explode the core of the star, accretions stops. 
This is the feedback component in the model. 

4.2 Numerical Setup 

We use the FLASH code to solve the non-relativistic hydro- 
dynamic equations on an adaptive-mesh refinement with 12 
levels. The simulation is done in a 2D cylindrical coordi- 
nates on a grid of size 1.5 x 10 9 cm in each direction. The 
equation of state (EOS) i s the Helmholtz equation of state 
dTimmes & SwestyL 12000b used in FLASH. This EOS in- 
cludes contributions from partial degenerates electrons and 
positrons, radiation and non degenerate ions. For the initial 
co nditions of the amb i ent ga s in the core we used the results 
of iLiebendorfer et aD d2Q05|) who made a ID simulations of 
the core bounce. We map their results from about 0.2 s after 
bounce for our simulations. For gravity we use FLASH'S 
multipole solver with only a monopole term (I = 0). In 
addition we take the newly born-NS to be a point mass 
M = IAMq in the center. The cooling function is taken 
to be 



= 5x l 3o£io. 

9 



-gxlO^X^Tfoergg-^- 1 , (8) 



here Tio is the temperature in 10 10 K, p is the density in 
grams, and X nuc is the mass fraction of free nucleons. 
The first term is the contribut ion of electron positron pair 
annihilation (Itoh et al. , 1996). The second is due to elec- 
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tron /positron capture on free nucleons (lOian & Woosley . 
Il996h , 

The jets are injected in pairs, alternating between the 
two sides. For each 0.1 s there is a 0.05 s jet-launching 
episode in the upper half plan, followed by a 0.05 s jet- 
launching episode in the lower half plan. The jets are in- 
jected in directions of 20° — 50° relative to the symmetry 
axes, and with a velocity of Vj = 10 10 cm s _1 . 

4.3 Results 

The first jet is able to penetrate through the collapsing enve- 
lope. Two shocks are created: a reverse shock where the jet 
is shocked, and a bow shock expanding outward shocking 
the core material. The bow shock accelerates the gas out- 
wards, while material more or less perpendicular to the jet 
initial direction continues to flow inwards. The two shocks 
heat the jet and the core material, respectively, to tempera- 
tures of ~ 2 — 7x1 9 K creating hot bubbles. In this range 
of temperatures the neutrino cooling is not significant and 
most of the energy is deposit inside the surrounding gas. Nu- 
clear reactions should occur inside the hot bubbles but this is 



not t aken into account in the simulation ( see lPapish & Soker 
2012l for discussion of r-process in these hot bubbles). 

The distance the first jet can penetrate depends on its 
properties, as well as on how fast it jitters. In our simulations 
this distance is about few x 1000 km. In the next episodes 
the jets are injected into a more dilute environment left by 
the expanding gas from previous episodes. This allows the 
jets to penetrate faster and catch up with previous shocks. 

This mechanism, were each time a jet in a different di- 
rection is injected, cause the hot bubbles to coalesce and 
form a more spherically shaped outflow; see lower left panel 
of Figure[l2l As the angles of the jets' jittering are changing 
within only 30°, the material in other directions keeps flow- 
ing inward to the NS. Accretion is reversed in all directions 
at ~ 3000 km, namely, where the hot bubbles coalesce into 
one large bubble that accelerates the core material outward. 
The evolution is seen in the first 3 panels of Figure fT2l 

There are some differences between the behavior of the 
jets in the supernova case and in the cluster case. These dif- 
ferences are mainly because the jets in supernovae are much 
denser than their surrounding material, while in the clusters 
the jets are less dense than their environment. 

The highly dense jets in our simulations of CCSN cause 
the outer froward shocks running to the core to be very fast. 
Hence, the forward shocks and the reverse shocks (where 
the jets are shocked) have similar velocities, implying sim- 
ilar post-shock temperatures as seen in Figures [T3l and ITU 
In clusters of galaxies and in planetary nebulae the jets are 
shocked to a temperature of two-three orders of magnitude 
higher than the temperature of the surrounding gas. 

Another difference is in the structure of the cocoon. A 
high density jet in CCSN flow creates a weak cocoon com- 
pared to the cocoon seen around jets in our simulations of 
clusters of galaxies. However, a large vortex exists inside the 



-10 -5 5 
1000 km 

Fig. 12 Comparing results from 2.5D simulations of jit- 
tering jets in core-collapse supernovae (Papish & Soker 
2012, in prep.) to an Ha i mage of the PN He 2-47, a low- 
excitation PN (image from lSahail20 00). Panels 1-3 show the 
log density (in cgs) of the inner most part of a 15M star 
after bounce. The times are after the first jet is injected. The 
similarity between the CCSN simulation and the PN image 
is in the structure where the bubbles protrude in different 
directions. 




-15-10 -5 5 10 15 
1000 km 



Fig. 13 Left panel: The temperature in the hot bubble at 
0.58 s after the start of the jets injection (Papish & Soker 
2012, in prep.); there are 6 episodes of jets in the upper side 
and 6 episodes of jets in the lower side. Right panel: HST 
image of the planetary n ebula Hu 2-1 obtaine d in the [N 
II] 65 8 3 line (image from lMiranda et al.1 ( 2001 )). The simi- 
larity between the CCSN simulation and the PN image is in 
the structure where bubbles are seen to catch up with previ- 
ously inflated bubbles. 
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Fig. 14 Upper panel: The temperature in the hot bubbles at 0.14 s after the start of the jets injection after one episode 
of jets in the lower side and two episodes of jets in the upper side. The contours represent the mass fraction of the jets 
material. Right upper panel: Zoom-in of the white square in the left panel. The velocity is shifted to the reference frame of 
the reverse shock. The velocity arrow in the white box is scaled to 10, 000 km s _1 . Lower panels: the results from Figure 
[7] panel 2, and Figure 0] panel (b) shown here again for comparison. 



bubble, similar to those in cluster simulations. To present 
the vortex, in Figure [141 we draw the velocity relative to the 
rest frame of the reverse shock. A clear counter-clock wise 
vortex is seen in the inset. The two lower panels show vor- 
tices from cluster simulations. The lower-left panel is the 
vortex of the flow shown in panel 2 in Figure [71 while the 
lower-right panel is the vortex of the flow shown in panel 
(b) in Figure HI 

Finally we note the very interesting similarity in mor- 
phologies of our simulated CCSN bubbles to those in some 
PNs, as shown in Figures [121 and \13\ The structure in our 
simulations of CCSNe is changing from a bipolar structure 
into a more spherical structure in about ~ 0.5 s. The PNs 
are relatively young. In old PNs such delicate bubble struc- 
tures are smeared out. 



5 JFM inside a common envelope 

Multidimensional numerical simulations of the common en- 
velope (CE) process usuall y include the gravity of the star 
and the companion (e.g.jLombardi et al.ll2006l: IPassv etal] 
1201 ll: iRicker & Taaml 120081 l20l5 1sandqui st et alJ Il998h . 
One of the main goals of these simulation s is to find the 
manner by which the envelope is ejected. ISokerl ( 2004b|) 
suggested that in some cases accretion disks are formed 
around WD and NS companions spiralling-in inside a CE. 
In these cases the accretion disk might launch two jets. Due 
to the very high velocity of jets from WDs and NS, the tem- 
perature of the post-shock jet material is much larger than 
that of the surroundi ng enve lope. Hot bubbles are formed. 
These bubbles. ISoke r ( 2004b) argued, can facilitate the ejec- 
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tion of the AGB envelope. We here conduct a preliminary 
study of this process for a WD orbiting inside the envelope 
of a massive AGB star. More realistic would be to study 
the jets blown by a NS, as it is not clear an accretion disks 
can be formed around s WD in a CE. Simulating jets from 
NS requires more sophisticated ingredients in the numerical 
code, and is left for a future study. 

The AGB star has a total mass of Magb = 5M , a core 
mass ofM core = O.77M , and a radius of Ragb = 3KLR . 
The density profile in the envelope is taken to be p ex r -2 . 
The relative velocity of the WD to the envelope is taken to 
be 50 km s _1 , somewhat lower the the Keplerian velocity 
due to envelope rotation. In this preliminary study we ne- 
glect the gravity of the WD, so that accretion is not treated 
here. This will have to be added in future simulations. The 
orbital separation is taken to be a = 150 Rq. 

The simulations are performed using the high-resolution 
multidimensional hydrodynamics code FLASH 4.0b, using 
a 3D uniform grid with Cartesian (x,y,z) geometry. We 
treat a small volume of the AGB envelope, and simulate the 
flow in a box with dimensions of 6Oi? x 30Rq x 30Rq 
in the x, y, z directions, with 512 x 256 x 256 cells along 
each axis (single cell size is~ O.IRq). We set the x axis 
along the relative velocity of the WD and the AGB en- 
velope, and inject the AGB envelope gas with velocity 
v x = 50 km s _1 from the minimum x boundary (x = 0) 
through the y — z plane. We impose outflow boundary con- 
ditions on all the other sides of the simulation box. The y 
axis is taken along the radial direction, such that the grav- 
ity is taken along the y axis and the density is p{y) — 
3.1 x lO 19 (?/+12Oi? )~ 2 g cm" 3 . The temperature is taken 
to be poly tropic with T ex p 2 / 3 , and with a value that en- 
sures hydrostatic equilibrium along the y axis. There is a 
mirror symmetry about the z = plane, and we apply a re- 
flective boundary condition there. A schematic drawing of 
the computational grid is presented in Figure [T5l 

We model the half space z > 0, and inject the jet in 
that side with a velocity of v z = 3000 km s _1 . The jet is 
launched from a square of 9 cells (3 x 3) in the plane z = 0. 
The mass loss rate of one jet is M\ = 2.5 x 1O _5 M yr _1 , 
such that the kinetic power of the two-jets is about the Ed- 
dington luminosity of a 1M© WD. Our main results are pre- 
sented in Figure [161 

In the present study we take the source of the jets to be a 
WD, but the results are relevant to accreting NS as well. As 
with the rest of the present study, we focus on morphologies 
of bubbles and the velocity structure in the bubbles and their 
surroundings. 

1 . Morphology. As expected, the AGB envelope bends the 
jet. This morphology is very similar to the one ob- 
served when radio jets are bent when a radio galaxy 
is moving through th e intra-galactic medium, e.g., 
O'Dea&Owenl (1986). These are called narrow-angle- 




tail radio sources. T he bending of radio jets w as simu- 
lated in the past, e.g jBalsara & Normanl (I1992|) . We ba- 



sically see a morphology similar to that of smoke flow- 



Fig. 15 A not-to-scale drawing of the WD motion within 
the AGB envelope. Our computational domain is marked by 
the box. The x direction is opposite to the velocity vector of 
the WD, v. The y direction is along the radial direction from 
the center of the AGB star. The fast jet is injected in the z 
direction. 



ing out of a chimney in a wind. However, we neglected 
convection in our simulation, something not completely 
justified as AGB envelopes possess a vigorous convec- 
tion. The typical size of the convective cells in AGB 
stars is the pressure scale height of l p (r) ~ 0.3r, which 
in our case is more or less the size of our grid. The con- 
vection velocity is close to the sound speed. With these 
typical values, we expect the convection to substantially 
influence the morphology. 

2. The flow near the jets' source. In addition to the 
shocked jets' material that is dragged by the envelope 
to large distances, some jets' material is flowing just 
behind the source. This flow can be seen in the panel 
(g) of Figure [T6l and w as studied in some detail by 
iBalsara & Normanl (1992). The hot and low-density re- 
gion formed by this flow can substantially affect the ac- 
cretion process onto the WD. 

3. Bubble fragmentation. As best seen in the two den- 
sity panels at t = 4.4 days, a separated bubble of 
hot shocked jet's gas is formed. This fragmentation is 
caused by instabilities, mainly the Kelvin-Helmholtz in- 
stability, as we already discussed for bubbles in clusters 
in section IX2l and can be seen in Figure HI 

4. Buoyancy. Hot bubbles are known to buoy outward in 
clusters of galaxies. The buoyancy of the hot bubble is 
seen in panel (f) of Figure [T51 as a motion in the +y di- 
rection. Gravity goes down in this panel, and the bubble 
buoys upward. 

5. Feedback: the AGB envelope. ISoked d2004bh discussed 
the influence of the hot bubbles formed by the jets on 
the AGB envelope. Although there is some energy in- 
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Fig. 16 The first three images (a,b,c) with the red background show the density maps at three different times given in 
days in the z — x plane passing through the WD source of the jet. The jet is injected from (x, z) = (10 12 cm, 0) in the 
+z direction. The envelope is injected from the left, in the +x direction. The density scale is in units g cm -3 . Panel (d) 
shows the temperature in the same plane in units of K. Panel (e) presents the velocity magnitude map in the same plane. 
Enlargement of the marked part is shown in panel (g). Panel (f) shows the density in the z = 0.5 x 10 12 cm plane. 
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put from the bubbles, jets from WDs are not expected to 

influence much the AGB envelope (maybe only in the 

very outer envelope regions). This is also seen in our 

simulation by the small volume of the hot bubbles. The 

volume of the bubble is smaller than the typical volume 

of a convective cell. However, NS can accrete at a very 

i - 1 1 — —i J 

high rate due to neutrino cooling (IChevarieru l993). and 

if they blow jets they can substantially disturb the enve- 
lope. 

6. Feedback: the accretion process. The flow behind the 
jets' source, as discussed in point (2) above, forms a hot 
region close to the source. This might substantially re- 
duce the accretion rate, hence forming a negative feed- 
back of the accretion process. This reminds us of the 
feedback in cooling flows. This feedback will be studied 
in the future when gravity is added to the jets' source. 

6 Summary 

A generic feature of many gravitationally bound systems 
is the tendency of compact regions to become smaller, and 
of the more extended regions to expand. The fundamental 
reason is that there is no real equilibrium state, and entropy 
(or another similar function, sometimes called H-function) 
can in principle increase unlimited. The gravitational energy 
released by a contracting small volume can substantially, 
and even catastrophically, affect the evolution of extended 
regions. 

There is the question of how the energy is transferred 
from the compact region to the extended regions. In stellar 
systems the energy transfer is done by gravity itself. Ex- 
amples include a triple system where the orbit of two stars 
shrinks and the tertiary body is ejected, and globular clusters 
where the core shrinks and the outer region expand. 

In most systems where the extended regions are filled 
with gas, radiation is inefficient in depositing energy to the 
gas. In cluster of galaxies and in galaxies, the hot gas is 
transparent to the radiation emitted by the AGN. As well, 
PNe absorbs a small fraction of the central stellar luminos- 
ity. In CCSNe at most ~ 1% of the energy carried by neu- 
trinos is absorbed. 

More efficient in influencing the ambient gas are jets 
launched by the compact object. The systems studied here, 
or to whom the results are relevant, and their order of mag- 
nitude properties are listed in Table [T] In some of these the 
interaction of the jets with the ambient gas is regulated by a 
feedback mechanisms (JFM). The emphasis of this paper is 
that despite the large differences between the systems listed 
in Table [TJ there are intriguing similarity in the interaction 
of jets with the ambient gas. These similarities can be used 
to learn from one system on another. For example, bubbles 
are observed directly only in CF clusters and in PNe. The 
direct observations can be used to compare the numerical 
simulations with observations. 

Our results strengthen the jet feedback mechanism 
(JFM) as a common process in many astrophysical objects. 



Part of the simulations were performed on the TAM- 
NUN HPC cluster at the Technion. The software used in 
this work was in part developed by the DOE NNSA-ASC 
OASCR Flash Center at the University of Chicago. This re- 
search was supported by the Asher Fund for Space Research 
and the E. and J. Bishop Research Fund at the Technion, and 
the Israel Science foundation. 
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